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Q Abstract 

^ I The impact of the presence of risk of destructive event on the silvicultural practice of a forest stand is 
d ] investigated. For that, we consider a model of population dynamics. This model has allowed us to make the 
S comparison without and with risk, and highlight the influence of the presence of risk of destructive event 
■ on optimal thinning and optimal rotation period. 
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^ ; Introduction 

^ I In terms of forest management, the first question that arises is : what is the optimal rotation period ? In 
the case where a calculation method to predict earnings for various rotation period is available, Faustmann 
^ (1849) proposed a formalism based on the expected discounted income. Many authors have successively 
'. improved or reformulated it ( [Ohlin ,1921], [Pearse, 1967], [Clark, 1976]). The risk of destruction has 
^ ' been introduced to forest stands by Martell (1980) and Routledge (1980) in discrete time. Thereafter, 
Reed (1984) has studied the optimal forest rotation in continuous time with the risk of fire. Thorsen and 
Helles (1998) analysed endogeneous risk. Buongiomo (2001) proposed a generalization of Faustmann 
approach using Markov Decision Process Models. Peyron and Heshmatol-Vezin (2003) were interested in 
natural risks incurred by forests in discrete time. More recently, Goodnow et al. (2008) took into account, 
in the stand management, that the thinning regime affect the proportion of standing trees damaged for 
ice-damaged, Amacher et al. (2008) highlighted the influence of silvicultural practices of a landowner on 
damage if fire occurs. 

Moreover, many authors have been studying ( [Naslund, 1969], [Schreuder, 1971], [Clark, 1976], [Kao 
and Brodie, 1980], [Roise, 1986], [Haight et al., 1992]), the determination of optimal thinning and cutting 
age, under certainty, using whole-stand models. More recently, Hyytiainen and Tahvonen (2003) studied, 
among other elements, the influence of the rate of interest on the rotation period, Cao et al. (2006) analysed 
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the effect of initial stand states on optimal thinning regime and rotation. Touza et al. (2008) investigated 
the ecology impact and interaction with the management. 

For the absence of risk of destructive events, all the production cycles are carried out to the same cutting 
age. When the risk of destructive event exists and is taken into account, we assume, as most authors cited 
above that the operator systematically decides to interrupt the current cycle and begins a new cycle. The first 
question that follows is about the impact of presence of risk of destructive event on silviculture. A second 
question, which is linked to the first one, is : what is the consequence of silviculture over the rotation 
period ? To take into account thinning in the calculation of the land value we used a model of population 
dynamics. Contrary to Xabadia and Goetz (2010) works based on age- structured models, we consider a 
simplified model of average individual type i.e. a model based on characteristics of medium trees, in order 
to facilitate interpretation and analyse more specifically the impact of a risk of destruction on silviculture. 
The study generalizes easily to more complex and realistic models. 

In a first part, we determine the land value without or with the presence of risk. In a second part we 
discuss the silviculture first in the reference limit case where individual tree growth is independent of tree- 
density (no density dependent growth) and then where individual tree-growth depends on the density (den- 
sity dependent growth). Under the assumption of no density dependent growth and for a fixed period of 
rotation, we compare the results in the absence and presence of risk. Then we deduce the impact of the 
presence of a risk to forestry. In a third part, under the assumption of density dependent growth and using 
the optimal rotation period, we simulate a Eucalyptus stand. This study, by taking into account explicitly 
implemented thinning, allows us to deepen the results obtained by Reed (1984) and clarify the underlying 
assumptions (not taking into account thinning and cleanup costs proportional to the damage). The interest 
to introduce a model of population dynamics and to use the proposed method is justified by the possibility 
of taking into account thinning and clearing costs depending on the severity of damage therefore. 

The land value 

In the first part we study the land value without the risk of destructive event through a model of popula- 
tion dynamics. Then we study the same land value with the risk of destructive event. 

Without the risk of destructive event 

We first consider a stand without the presence of risk of destructive event. The study of this case will 
allow us to define a benchmark management of the stand. 

For a cutting age, or rotation period T and a rate of thinning by unit of time /i(.), the land value Wq (up 
to a constant ci) is the discounted value of cutting incomes net cost of replanting : 

1=1 

where V(/i(.),T) is the income generated by the cutting at time T and ci is the cost of replanting. 
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V{h{.), T) is by definition the sum of the thinning income on period [0, T] and the income at final cutting 
age T. 

The land value can also be interpreted as, Wq is the instantaneous value of income from the forest in 
time T discounted at the initial time and is solution of : 



Wo = {Wo + V{h{.),T)-c,)e-'^ 

A model of population dynamic 

To take into account the thinning in the calculation of the land value, i.e. to express V(/i(.),T), we 
introduce a model of population dynamics. The considered model is an average tree model : the state 
variables are the number n of trees per hectare and the averaged tree-basal area s measured at breast height 
(1.30 meter from the ground). The evolution of these two variables is govemed by the system of ordinary 
differential equations : 



^ = -(m(t) + h{t))n{t) 
ds{t) 



dt 



G{n{t),s{t),t) 



where m(.) is the natural mortality and G{., ., .) is the possibly density dependent growth function : 
individual tree-growth depends on the tree density n{t). 

The considered population dynamic model permits to explicitly manage the trees. The use of more 
complex and realistic models with different classes of tree basal areas, instead of the model presented, 
would not cause any additional methodological problem. In order to facilitate interpretation and analysis of 
properties and results, we chose to assume identical the characteristics of the trees. 

Link with stand models 

Such models can be linked with stand models in the following way, let the basal area of the stand 
S = n.s and assume that nG{n, s,t) = g{n, s)T{t) and that g{n, s) = go{S) depends only on 5" = n.s. 
Then the dynamics of the basal area is given by : 

^ = go{S{t))m - (m(t) + h{t))S{t) 

In the particular case of m = and of the price p proportional to s, by substituting first the basal area to 
the volume and secondly hS to h we obtain the Clark's model ( Clark, 1976 pp 263-269). 

The forest income 
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Once chosen the model of population dynamics we can express the total income. Total income V{h{.),T) 
for fixed thinning h(.) is the sum of the thinning income H(n(.), h{.),T) on the period [0, T] and the 
final income Vo{n{T), s(T)) : 

V{h{.),T) = H{n{.), s{.),h{.),T) + Vo{n{T),s{T)) 
The thinning income H{n{.), s{.), h{.), t) on [0, t] actualized to time t is : 

H{n{.),s{.),h{.),t) = [ j9(s(M))/i(w)n(w)e*(*-")(iM 



^0 

where n(.) and s(.) are solutions of the dynamic model. 
Calculation of the land value 

The land value is given by : 



V{h{.),T)-c^ 
- e^^ - 1 

fT 

with : V{h{.),T) = I p{s{u))h{u)n{u)e^^'^-''Uu + Vo{n{T), s{T)) 



The maximal value of the land value defined by Faustmann is obtained by solving the problem : 

V{h{.),T)-c, 



(Vfo) ■ max 



h{.),T e^^ - 1 

The maximization of the Faustmann value Wo taking into account thinning h{.) for the cutting age T 
can be decomposed in two steps : first we maximize V{h{.),T) with respect to h{.) : 

{So) Vo(T) =maxV(/i(.),^) 

hi.) 



then we maximize ° ^ with respect to T with Vo(T) resulting from the first step. 



In the presence of risk of destructive event 

Following Reed (1976), we suppose that destructive events occur in a Poisson process i.e. that destruc- 
tive events occur independently of one another, and ramdomly in time. 
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The distribution of the destructive event time is an exponential with mean — : F{x) = 1 - e-^"" where A 

A 

is the expected number of destructive events per unit time. No assumption is made on the type of destructive 
events. Destructive event can also be a storm, a fire, a disease or an insect infestation. We assume that the 
event is always of the same type. 

We assume that 9t is the proportion of damaged trees following a destructive event and 9' is the rate 
depreciation of timber due to an influx of wood on the market. 9 and 9' are random variables positively 
correlated : if 6' = (no damage) then 9' = and more 9 is high, more 9' is also. Depending on the type 
of risk, 9t and thus 9[ may be conditionally distributed to n{.), s{.), to thinning h{.) or to another control 
parameter. We define the expectations a of (1 — 9^) and ap of (1 — — 9[), due to the assumptions 
on the random variables 9t and 9'^ we deduce ap < a. In order to lighten the presentation we omit the 
n{.), , h{.) dependancy of a and ap. The modelling is suitable for different types of destructive risk : the 
specificity is reflected in the choice of random variables 9, 9' distribution laws and therefore of functions a 
and ap. 

For a cutting age T and a fixed thinning h{.), the land value is the actualized value at initial time of 
the sum of two terms. The first one is the expectation (with respect to the time of event) of the sum of the 
land value and the expectation (with respect to 9t and 9[) of the total income Vi(/i(.), t) for the period [0, t] 
minus clearing costs C{n(t) , t) in case of a destructive event at time t. The second one is the sum of the land 
value and the total income V(/i(.), T) for the period [0, T] in case of no destructive event. Then is the 
solution of the following equation : 



Wo= / [iyo + Vi(/i(.),t)-Ci-C(n(t),t)]e-^*rfF(t) + (iyo + V(M.),3^)-Ci)e-^^(l-F(r)) (1) 



where Vi(/i(.),t) = E{V{h{.),9,9\t)) with V(/i(.), ^, t) the total income for ^ and C(n(f), t) = 
E{Cn{9t, n{t))) with Cn{9t, n{t)) the clearing costs for 9t. 

The forest income 

V(/i(.), 9, 9' , t) is the sum of the thinning income H{n{.), s(.), h{.),t) during [0, t] actualized at time t 
and the total income Vf{9, 9' , t) at time t : 



The final income VF{9,9',t) is assumed proportional to the final income without risk of destructive 
event K)(n(r), s(T)) and is given by : Vf(^, 9\ t) = {1- 9){1 - 9')Voin{t), s{t)). From definition of a and 
ttp we deduce the total income expectation : 




V(/i(.), 9, 9', t) = H{n{.), s{.),h{.),t) + Vf{9, 9', t) 



V,{h{.),t) = E{V{h{.),9,9' ,t)) = H{n{.),s 



{.),h{.),t) + ap{t)Vo{n{t),s{t)) 
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In Reed (1976), the total income generated in case of risk of destructive event is assumed proportional 
to the income generated in the case without risk, which is not longer true if one wants to take the thinning 
into account. 

The clearing costs 

The clearing costs C„ are of two types : the first one concerns the damaged trees 6tn{t) and the second 
one concems the survival trees (1 — 9t)n(t). The respective contribution of each type of cost depends on 
the type of destructive event. Assuming a linear dependency ( Amacher et al., 2008), we then deduce the 
expression of the clearing costs : 

CniOu nit)) = C2 + cAnit) + c,(l - ^t)n(t) 
From the definition of a, we deduce the clearing costs expectation : 

C{n{t),t) = E{Cn{et, n{t))) = C2 + Cn{t)n{t) 

where c.„ is defined by c„(t) = Crf(l — a{t)) + Csa{t). Thus, min(Q, Cg) < Cn{t) < max(Q, Cg). 

Except for the fire risk ( Amacher et al., 1976), the clearing costs mainly concems the damaged trees. 

Remark : Reed (1976) assumed that clearing costs were independent of the number of damaged trees. 
Using a model of population dynamics allows us to take into account the clearing costs proportional to the 
damage. 

The salvageable function a 

The salvageable function a dependency is specific of the type of risk of destructive event. Hence : 

- for a storm risk, following Hanewinkel (2008), a is a decreasing function of the ratio height-diameter 
^ (where d is the tree-diameter d = \J^) and the time t : then a can be written as a{-^^=, t). More 
generally a is a function of n(t) and s(t). 

- for ice damage risk, following Goodnow et al. (2008), a depends on the planting density uq and on 
the rate of thinning h{.). 

- for a fire risk, following Amacher et al. (2008), a depends on the planting density uq and on the level 
of preventive intermediate treatment effort z{.). In that case, 2;(.) is an additional control parameter. 

Calculation of the land value 

From (1) and expression of Vi(/i(.), t) we deduce the land value : 
Proposition 1 : In the presence of the risk of destructive event the land value is given by : 
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6 + X V,{h{.),T)-c, A 

w/zere ; Vi(/i(.), = / b(s(t))/i(t)n(t)+Aap(t)yo(n(t), s(t))-Ac„(t)n(t)]e(^+^)(^-*)dt+\/o(n(r), s(T)). 
Proof: From (1) and expressions of V(/i(.),t) and Vi(/i(.), t) : 



Wo = + H{n{.),s{.), h{.),t) + ap{t)Vo{n{t),s{t)) - c, - C2 - Cn{t)n{t)]e~''dF{t) 

Jo 

+ {Wo + H{n{.), s{.),h{.), T) + Vo{n{T),s{T)) - ci)e-^^(l - F{T)) 
We then deduce the land value : 

... ^ + A Jo(n(.), h{.),T) + HKr), g(T)) - ci A 

^0 = ^ e(^+A)T _ 1 ^(Cl + C2) 

where Jo (n (.), s (.),/;.(.), T) has the following value : 

^[H{n{.), s{.),h{.),t) + a,{t)Vo{n{t), s{t)) - c^{t)n{t)]e'^^-'^de'^^~'^ + H{n{.), s{.),h{.), T). 

After changing the integration order in / H{n{.), s{.), h{.),t)e^^'^~^^de^^'^~*\ Jq becomes : 

Jo 

Mn{.),s{.), h{.),T) = f [p{s{t))h{t)n{t) + \a^{t)Vo{n{t) , s{t)) - Ac„(t)n(t)]e(^+^)(^-*)6?t 

then the result. □ 







As in the case without risk, we find that the land value can be deduced from the income Vi(/i(.), T). The 
following differences occur : first the rate of interest 5 is replaced by 5 + A, secondly V(/i(.), T) is replaced 
by a modified expression of the income Vi(/i(.), T). The second difference is reflected by the substitution 
in the case without risk of term H{n{.), s{.),h{.),T) by the term Jo(n(.), s(.), T), or even more 
precisely, the substitution of by p{s{t))h{t)n{t) + X{ap{t)Vo{n{t) , s{t)) — Cn{t)n{t)). 

The maximal value of the land value is obtained by solving : 

: max Wq = — ^— 7 Ci + C2 

h{.),T 5 e(*+A'^ - 1 6 
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As in the case without risk, the maximal value of the land value with respect to the thinning h{.) and the 
cutting age T can be decomposed in two steps : first we maximize Vi(/i(.), T) with respect to h{.) : 

(Si) Vi(T) = maxVi(M.),^) 

h(.) 

then we maximize ^ ~^ — ~~ — fI"^! + ^2) with respect to T with Vi(T) resulting from the first 
step. 

Silviculture for a fixed cutting age T and a no density dependent growth 

We consider first the limiting case of no density dependent growth which will be used as a reference in 
the case study of density dependent growth. 

Without risk of destructive event 

Let us consider the case where individual growth is not density dependent. In this case the evolution of 
the tree-basal area s does not depend on the tree-number n, then does not depend on the silviculture. 

We study the maximization of the land value with respect to thinning h{.) for a fixed cutting age T : 

{So) maxV(/i(.),T) = /" p{s{t))h{t)n{t)e^^^-'^dt + Vo{n{T), s{T)) 
hi-) Jo 

with the constraint < h{t) < h. 

Assume that the final income is given by : Vo{n, s) = p{s)n. From the no density dependence of 
the individual growth, p{s{t)) is independent of the thinnings and only depends on t. Then we define 
R(t) = p{s{t)). We denote the functions ttq and Ho : 7ro(t) = R'{t) - {6 + m{t))R{t) and Ilo{t) = 

e^^ ('^+™("')+'^)'^"'7ro(n)(iu. Applying the maximum Pontryagin Principle to the problem Vq (see Appen- 
dix A) we can deduce the proposition : 

Proposition 2 : Assume ttq is decreasing, the individual tree-growth is not density dependent. Consider a 
fixed cutting age T, then the optimal thinnings are given by : 

- ifnoiT) > then K = 

- if TXo{T) < then it exists < < T such that h^,{t) = for t < t^, and h{t) = h pour t > t^,. 
Moreover ifUo^O) > then = 0, else is the unique solution of: 
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Remark : The commutation time is a function of the cutting age T. 
In the presence of risk of destructive event 

We study the maximization of the land value with respect to thinning h{.) for a fixed cutting age T : 



with the constraint < h < h. 

In this paragraph, we consider a function a depending only of time t. The chosen function a can be 
interpreted as the expectation with respect to the state variables of a more precise function depending on 
the state variables. 

In case of a storm risk, a is a function of the ratio height-diameter and time t. For a no density dependent 
growth, the tree-basal area s and hence the diameter is an explicit function of time t, then a depends only 
on t. In this case, the approximate function is exact. 

Due to the difference in the criterion without and with risk, we deduce that, for the fixed rotation period 
T, the silviculture differs and depends explicitly on A, a and ap. So we will pay attention to the consequence 
for the silvicultural practice. 

We consider, as in the case without risk, a no density dependent growth for the trees to facilitate the 
comparaison. 

Let denote the functions ttx and 11^ : T^xit) = R'{t) — (A(l — ap{t)) + 6 + 'm{t))R{t) — Ac„(t)) and 

^x{t) = e-^" (^+'^+™("')+'^)'^"'7rA(ti)dM. Applying the maximum Pontryagin Principle to the problem V\ 
(see appendix A) we deduce the proposition : 

Proposition 3 : Assume tx\ is decreasing, the individual tree-growth is not density dependent. Consider a 
fixed cutting age T, then the optimal thinnings are given by : 

- ifixxiT) > then /i* = 

- if 7rx{T) < then it exists < < T such that h^{t) = for t < and h{t) = h for t > t^,. 
Moreover ifllx{0) > then = 0, else is the unique solution of: 



Comparaison : without and with presence of risk 

By comparing the without and with risk criteria, we remark that the rate of interest 5 is replaced by 5 + A 



maxVi(M.),T) = Jo(n(.),s 



{.),h{.),T) + Vo{n{T),s 
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and the rate of thinning h{t) is replaced by h{t) + X{ap{t) " ) if Vo is defined by Vo{n, s) = p{s)n. 

p[s[t)) 

Moreover, from the definition of tta, we deduce that tta is a decreasing function of A. Thus, at least, in 
the vicinity of the cutting age T, the greater A, the more frequently h(T) will be equal to h in the following 
sense : for Ai < A2, with the associated thinnings hi and /12, if hi{T) = h then h2{T) = h. 

By comparing the two propositions, without considering the clearing costs, the natural mortality m{t) in 
the case without risk is replaced by the mortality due to events m{t) + \(1 — ap{t)) in the presence of risk. 
It is equivalent also, from a mathematical point of view, to replace the fixed discount rate 5 by the variable 
discount rate 5 + A(l — a!p(t)) in the previous problem. 

By comparing the results of the two propositions we deduce that, for a fixed rotation period T, it is 
usually best to do thinning at least at the end of the period in the presence of risk. This comparison confirms 
the previous property deduced from the definition of tta. 

Comparing the results of the two proposals is permitted if the rotation periods are identical. If we 
consider the maximization problem, with respect to the rotation period, the rotation periods have no reason 
to be the same. In that case the comparison is not permitted and only simulations can allow us to compare 
the respective thinning. We will therefore perform simulations. 

Silviculture for a density dependent growth 

We now consider the case where individual growth is density dependent. If the growth is weakly density 
dependent, by continuity with the case of no density dependent growth, the obtained results are still valid at 
fixed rotation period T. For a greater density dependence, if {G{n, s, t)ny^ is sufficently small (see Appen- 
dix B), the optimal thinnings are the same as in the previous case in the vicinity of T. If this is not the case, 
we cannot obtain analytical results for the solutions, then simulations are required. 

We are interested in a stand of Eucalyptus. In the absence of more precise information on the structure 
of the salvageable function a for a Eucalyptus stand, we restricted our analysis to a constant. The function 
of individual growth is given by : 

, 0.7445(1 - e-o-482nsN ^^/^n 

G{n,s,t) = —— 

n dt 

^ 

where H{t) is the high at time t : H{t) = Hq{1 — e and Hq = 30 the limited high which depends 
mainly on soil fertility ( Saint Andre et al., 2002). The structure of the growth function is generic and can 
be used for other species. 

For the clearing costs, we neglected the clearing costs of the second type (c„(t) = q(1 — a{t)). 

The weight of the trees of basal area s and high h is given by : v{s, H, t) = 0.29+(127.8+0.32t)si/inkg 
( Saint Andre et al., 2005). The price is assumed to depend on the weight : p{s, t) = O.lf (s, H{t),t) — 0.25. 
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The determination of the cutting age T is important because of its impact on silvicuhure. To better 
describe the silviculture in the presence of random risk, it is wiser to look at the effective cutting age T and 
the effective final tree-basal area S. Thus we calculate the respective expectations and variances : 



E{r) = [ tdF{t) + T(l - F(T)) = 



^0 





AT) + 



F2 



(T) 



■T 



■T 



E{S) = / s{t)dF{t) + s(T)(l - F{T)) and Var{S) = / s\t)dF{t) + s\T){l - F{T)) can be 



derived from the simulations. 

Results and Discussion 

The unit of time for the rotation period is the month. We suppose : m = 0.0042 month^^, A = 0.0075 
month-\ S = 0.0034 month-\7i = 0.075 month-^ 

Assume first, that in case of destructive event, the destruction of stand is total (a = 0,ap = 0). With the 
constraint of no thinning, we found (Table 1) the classical well established result : risk implies a shortening 
of the optimum rotation period. By considering the optimal thinning, the optimal rotation period is larger 
than in the case without risk but the expected effective rotation period is smaller and remains of the same 
order of magnitude as in the case without risk. This can be explained by the fact that, in the presence of 
risk, the optimal rotation period T is achieved with a relatively low probability : 1 — F{T) = e~^^. To 
complete the study, the standard deviation of the effective rotation period was calculated. Its value doesn't 
vary and is about 20 months. The expected effective final tree-basal area varies only slightly depending on 
the scenarios. 

If the destruction is only partial i.e. a portion of the stumpage is salvageable (a = 0.6, ap = 0.4), 
without risk and for the two considered tree-density (n = 650 or 1650 stems/ha), there is no thinning for 
the optimal solution (Table 2). In presence of risk, with the constraint of no thinning, the optimal rotation 
period is close to the previous one. In contrast, if we optimize allowing thinning, the optimal rotation period 
is greater and thinnings are to be done. Similarly to the previous case of total destruction, the expected 
effective rotation period remains of the same magnitude as in the case without risk. 

The initial tree density in the studied range did not influence the qualitative behavior of optimal ma- 
nagement. Comparing Table 1 and Table 2 for a density of 650 stems/ha, we remark that, the greater the 
potential damage, the greater the cutting age will be and the earlier the beginning of the preventive thinning. 
From Table 2, we deduced that the greater the initial tree-density, the later the beginning of the preventive 
thinning will be and the lower the difference between the land value in presence of risk without and with 
thinning. In Table 3, for the commonly used cutting age value T = 84 months, we found similar properties 
for the different optimizations. 
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The presence of risk of partial or total destruction involves earlier thinning. Because of early thinning 
the amount of standing trees is smaller at time T. Thus, because of thinning, the rotation period T, can 
be extended and greater than the rotation period without risk. The earlier thinning provides a kind of self- 
insurance against risk. We therefore make endogenous the risk through optimization. 

Hyytiainen and Tahvonen (2003) showed, in another context under certainty, that a greater rate of in- 
terest may lengthen the optimal rotation for non optimal initial density. The density of 650 stem/ha in not 
optimal. We have showed that the risk involves in particular the substitution of the rate of interest by the rate 
of interest plus the expected number of events A, the consequence being an increase of the rate of interest. 
Thus our result is consistent with the results of Hyytiainen and Tahvonen. By observing the curves of the 
land value Wq (for a density of 650 stems/ha and a = 0.6, ctp = 0.4), depending on the rotation problem T 
with optimal thinning, we find that the land value least varies in the vicinity of the optimal rotation period 
with risk (Figure 2) than without risk (Figure 1). This is another consequence of the fact that, with risk the 
optimal rotation period T is achieved with a relatively low probability. 

We consider the case where the rotation period is determined by other considerations. We take the 
commonly used value T = 84 months (Table 3). Without and with risk, it is optimal to practice thinning. 
However, in the presence of risk, optimal thinning starts earlier. 

The results obtained in Reed (1976) are valid only under the following assumptions : the manager does 
not practice thinning and the clearing costs in case of destructive event are fixed. The possibility to consider 
thinning and clearing costs depends on dammage severity and therefore justifies the interest to introduce a 
model of population dynamics. 

Conclusion 

We have studied the management of a stand in the presence of risk of destructive event. In order to de- 
termine the optimal thinning relative to the Faustmann criterion, we have considered a model of population 
dynamics, the choosen model is of average tree type. This model has allowed us to make the comparison 
without and with risk and highlighted the influence of the presence of risk of destructive event on optimal 
thinning. 

Specifically, the obtained land values, without or with the risk of destructive event, highlighted diffe- 
rences in the criteria to be maximized. In the case of no density dependent individual growth, we have 
highlighted the impact of the presence of risk on the strategies, generically regardless of the considered 
species. 

In the case of density dependent growth, the calculations for a stand of Eucalyptus have shown that the 
presence of risk of destruction event involves earlier thinning and a greater rotation period, for the optimal 
strategy. 

The obtained results are conditioned by the choice of an individual tree growth model and by the speci- 
fication of a weight model and a price model of trees. Other studies using models adapted for other species 
would make the obtained results more generic. 



12 



Acknowledgements 

The author would like to thank 0111 Tahvonen for helpful discussions of the subject matter. The author is 
grateful to the anonymous reviewers, whose comments and suggestions helped improve the quality of this 
paper. 

References 

Amacher et al., 2008 Amacher, G.S., Malik, A.S. and Haight, R.G., Forest landowner decisions and the 
value of information under fire risk. Can. J. For. Res. 35 (2008), pp 2603-2615. 

Buongiomo, 2001 Buongiomo, J., Generalization of Faustmann's Formula for Stochastic Forest Growth 
and Prices with Markov Decision Process Models. Forest Science, 47(4) (2001), pp 466-474. 

Cao et al., 2006 Cao, T., Hyytiainen, K., Tahvonen, O. and Valsta, L., Effects of initial stand states on 
optimal thinning regime and rotation of Picea abies stands. Scandinavian Journal of Forest Research, 21(5) 
(2006), pp 388-398. 

Clark, 1976 Clark, C.W., Mathematical Bioeconomics, Wiley, New York (1976). 

Faustmann, 1849 Faustmann, M., Berechnung des Wertes welchen Waldboden sowie noch nicht haubare 
Holzbestande fiir die Weldwirtschaft besitzen. Allgemeine Forst-und Jagd-Zeitung, 25 (1849), pp 441-445. 

Goodnow at al., 2008 Goodnow, R., Sullivan, J. and Amacher, G.S., Ice damage and forest stand manage- 
ment. Journal of Forest Economics, 14(4) (2008), pp 268-288. 

Haight et al., 1992 Haight, R.G., Monserud, R.A. and Chew, J.D., Optimal Harvesting with Stand Density 
Targets : Managing Rocky Mountain Conifer Stands for Multiple Forest Outputs. Forest Science, 38(3) 
(1992), pp. 554-574. 

Hanewinkel, 2008 Hanewinkel, M., Storm Damage Modelling in Southwest Germany based on National 
Forest Inventory Data. The final seminar of the Stormrisk project. Boras, Sweden, 2008-10-02. 

Hyytiainen and Tahvonen, 2003 Hyytiainen, K. and Tahvonen, O., Maximum Sustained Yield, Forest Rent 
or Faustmann : Does it Really Matter! Scandinavian Journal of Forest Research, 18(5) (2003), pp 457-469. 

KaoBrodie Kao, C. and Brodie, J.D., Simultaneous Optimization of Thinnings and Rotation with Conti- 
nuous Stocking and Entry Intervals. Forest Science, Monograph 22 (Supplement to Number 3) (1980), pp. 
338-346. 

Martell, 1980 Martell, D.L., The optimal rotation of a flammable forest stand, Canadian Journal of Forest 
Research, 10(1) (1980), pp. 30-34. 

Naslund, 1969 Naslund, B., Optimal Rotation and Thinning. Forest Science, 15(4) (1969), pp. 446-451 

Ohlin, 1921 Ohlin, B.,) Concerning the question of the rotation period in forestry. Journal of Forest Eco- 
nomics, vol l,n°l-1995, (1921), pp. 89-114. 

Pearse, 1967 Pearse, PH., The optimal forest rotation. Forest Chron., 43 (1967), pp 178-195. 



13 



Peyron and Heshmatol-Vazin, 2003 Peyron, J.L. and Heshmatol-Vazin, M., La modelisation de la foret 
landaise : portee et limites. Communication 7th Conference ARBORA, Pessac, 14-15 Dec 2003, 17p. 

Reed, 1984 Reed, W.J., The Effects of the Risk of Fire on the Optimal rotation of a Forest. JEEM, 11 
(1984), pp 180-190. 

Roise, 1986 Roise, J. P., A Nonlinear Programming Approach to Stand Optimization. Forest Science, 32(3) 
(1986), pp. 735-748. 

Routledge, 1980 Routledge, R.D., Effect of potential catastrophic mortality and other unpredictable events 
on optimal forest rotation policy. Forest Science, 26 (1980), pp 386-399. 

Saint Andre, 2002 Saint-Andre, L., Laclau, J. P., Bouillet, J.P, Deleporte, P., Mabiala, A., Ognouabi, N., 
Bailleres, H. and Nouvellon, Y., Integrative modelling approach to assess the sustainability of the Eucalyp- 
tus plantations in Congo. Fourth workshop, lUFRO Working Party S5.01.04 « Connection between Forest 
Resources and Wood Quality : Modelling Approaches and Simulation Software*. Harrison Hot Springs 
Resort, British Columbia, Canada, September 8-15, 2002. 

Saint Andre, 2005 Saint-Andre, L., M'Bou, A.T., Mabiala, A., Mouvondy, W., Jourdan, C, Roupsard, O., 
Deleporte, P., Hamel, O., and Nouvellon, Y., Age-related equations for above- and below-ground biomass 
of a Eucalyptus hybrid in Congo. Forest Ecology and Management, 205 (2005), pp 199-214. 

Schreuder, 1971 Schreuder, G.F, The Simultaneous Determination of Optimal Thinning Schedule and 
Rotation for an Even-Aged Forest. Forest Science, 17(3) (1971), pp. 333-339. 

Thorsen and Helles, 1998 Thorsen, B.J. and Helles, F, Optimal stand management with endogenous risk 
of sudden destruction. Forest Ecology and Management, 108(3) (1998), pp. 287-299. 

Touza et al., 2008 Touza, J., Termansen, M. and Perrings, C, A Bioeconomic Approach to the Faustmann- 
Hartman Model : Ecological Interactions in Managed Forest. Natural Resource Modeling, 21(4) (2008), pp 
551-581. 

Xabadia and Goetz, 2010 Xabadia, A. and Goetz, R.U., The optimal selective logging time and the Faust- 
mann formula. Journal of Forest Economics, 16 (2010), pp 63-82. 

Appendix A. Solving the optimisation problem {Vx) with a no density dependent 
growth 

We consider, for a fixed rotation period T, the following problem {V\) : 

max / [R{t)h{t) + X{ap{t)R{t) - Cn{t))]n{t)e^^-^^^^^-'Ut + R{T)n{T) 
K) Jo 

with the constraint < h{t) < h. 
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We apply the maximum Pontryagin Principle, the Hamiltonian is 



H = [R{t)h + \{ap{t)R{t) - c„(t))]ne(''+^)(^-*^ - /i(m(t) + h)n 
The first order conditions are : 



h*(t) maximizes [i?(t)e(^+^)(^~*) - ii{t)]hn with h e [0, h] 

^ = = - mmt) + XMt)Rit) - c„(t))]e(^+^)(^~*) + fi{t){m{t) + h{t)) 

at on 

with the transversality condition fi{T) = R(T) 

We consider the evolution of the function l{t) defined by : l(t) = i?(t)e(^+^^(^"*^ - fi{t). 
dl{t) 



dt 



{R'{t) -{6 + A(l - apit)) + m{t))R{t) - Ac„(t))e(^+^)(^-*) + /(t)(m(t) + hit)) 



and define the function n and IT by : 7r(t) = R'{t) - [6 + A(l - ap{t)) + m{t))R{t) - Xcn{t) and 

We then deduce the following Proposition : 

Proposition A.l : Assume tt is decreasing, the individual tree-growth is not density dependent. Consider a 
fixed cutting age T, then the optimal thinnings are given by : 

- ifTr{T) > then K = 

- ifniT) < titen il existe < t^ <T such that h^{t) = Ofor t < t^^ and h{t) = hfor t > t^t. Moreover 
2/11(0) > then = 0, else t^ is the unique solution of: 

g/J{5+A+m(«')+/^)d«V(u)dM = 

Proof: From 1{T) = and considering the equation in backward time t' = T — t, we deduce that in the 

vicinity of T, and 7v{t) are of opposit sign, then —^r' ^^'^ '^(^) have the same sign : hence l{t) < 

(resp. > 0) if 7r(t) > (resp. < 0). If n{T^) > 0, 7r(t) > forall t then l{t) < forall t which implies 

h{t) = 0. If 7r(r,) < 0, from l{t) = -e- fu(^^''^+^)<i^u{t) we deduce that either n(0) < 0, l{t) < and 
l(t) cannot change of sign or 11(0) > and l(t) changes of sign for a unique time t*. 
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Appendix B. Solving the optimisation problem (Vx) with a density dependent growth 

We consider, for a fixed rotation period T, the following problem (Vx) : 

max r [p{s{t))h{t) + \{ap{t)p{s{t)) - c„(t))]n(t)e(^+^)(^-*»dt + p{s{T))n{T) 
K) Jo 

with the constraint < h{t) < h. 

We apply the maximum Pontryagin Principle, the Hamiltonian is : 

H = [p{s)h + X{ap{t)p{s) - c„(t))]ne(''+^)(^"*) - ^„(m(t) + h)n + fisG{n, s) 
The first order conditions are : 



h*{t) maximizes [p{s{t))e^^+^^'-^-^^ - ^{t)]hn with h G [0, h] 

= = - Witmt) + Ka,{t)p{s{t)) - c.(t))]e(^+^)(^-*)) + iin{t){m{t) + h{t)) 

- ^is{t)Gn{n{t),s{t),t) 

= = - {p'm)h{t) + A«,(ty(5(t)))n(t)e(^+^)(^-*) - fiMGs{n{t),s{t),t) 

with the transversality conditions /x„(T) = p{s{T)) and /is(T) = p'{s{T))n{T). 

We consider the evolution of the function l{t) defined by : l(t) = p(s(t))e^'^+'^^^^~*^^ - /i„(i). 



^ =(p'(s(t))G(n(t), .(t)) -{6 + A(l - a,(t)) + m(t))p(s(t)) - Ac„(t))e(^+^)(^-*) 

which allow us to give the following result : 

Proposition A.2 : Ifp{s) = Cs"- and a{G{n, s, t)n)'^ < {5 + A(l — ap{T)) + m(T))s(0) then the thinning 
h{t) = h is optimal in the vicinity of the cutting age T. 

dl 

Proof: From 1{T) = we deduce : — (T) = -{6 + A(l - ap{T)) + m{T))p{s{T)) - \Cn{T) + 
p'{s{T)){G{n, s, t)n)'^iT). From given conditions, -^(T) < then l{t) > in the vicinity of T. 
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Table 1 : Optimized land value with respect to T. a = 0, Op = 0. 



Scenario Optimal thinnings 

(month" ^) 

650 stems/ha 
Without risk 

ma.Xh(,)^T h = 

With risk 

no thinning 

maxy, no thinning 

max/i(,)^r h(t) = h,t> 36.5 



Cutting age Land value Expected effective 

cutting age 



(month) (euro) (month) 

58.5 2137.5 58.5 

58.5 658.1 47.3 

54 673.3 44.4 

69.5 810.4 54.2 



Table 2 : Optimized land value with respect to T. a = 0.6, ap = 0.4. 



Scenario Optimal thinnings 

(month" ^) 

650 stems/ha 
Without risk 

ma.Xh(,)^T h = 

With risk 

no thinning 

maxT, no thinning 

max/i(.)^r h{t) = h,t> 43.5 

1650 stems/ha 
Without risk 

max/i(.)^T h = 

With risk 

no thinning 

max^, no thinning 

max/i( h{t) = h,t> 46.5 



Cutting age Land value Expected effective 

cutting age 
(month) (euro) (month) 



58.5 


2137.5 


58.5 


58.5 


1108.9 


47.3 


57.5 


1109.8 


46.7 


65.5 


1149.0 


51.7 


59.5 


2497.2 


59.5 


59.5 


1230.7 


48.0 


58.5 


1232.1 


47.3 


64.5 


1251.1 


51.1 
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Table 3 : Land value for T = 84 months. 



Scenario Optimal thinnings Cutting age Land value Expected effective 

cutting age 

(month"^) (month) (euro) (month) 

650 stems/ha 
Without risk 

max/i(,) h{t) = h,t> 60 

With risk 

h{t) = h,t> 60 

max/i(.) h(t),t> 43.5 



84 

84 
84 



1914.1 

1018.4 
1099.3 



84 

62.3 
62.3 



h{t) = h,t> 64.5 



1650 stems/ha 
Without risk 

With risk 

h{t) = h,t> 64.5 

max/i(.) h{t) = h,t> 48.5 



84 

84 

84 



2224.0 

1099.3 
1194.9 



84 

62.3 
62.3 
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Fig 1. Without risl( case. Ttie abscissa sliows tlie rotation period T in monttis, ttie ordinate the land value 
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